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Abstract 

Monolayer transition metal dichalcogenides offer the possibility of optical control of the valley 
degree of freedom. In order to asses the potential of these materials in applications, detailed 
knowledge of the valley dynamics is essential. In this work, we apply low temperature time-resolved 
photoluminescence (PL) measurements to investigate exciton valley relaxation dynamics and, in 
particular, its behavior under strong excitation. At the lowest excitation powers the inter valley 
scattering time is ~ 50 ps, but shortens by more than a factor of two at the highest powers. 
We attribute this acceleration to either heating of the exciton system or the presence of a dense 
exciton gas, which could influence the exciton valley properties. Furthermore, we analyze the PL 
dynamics of excitons and trions. We find that the PL decays for all peaks are bi-exponential and 
approximately independent of the excitation power. We attribute the short decay to radiative 
recombination and escape to a reservoir of dark states. The long decay is ascribed to a transfer of 
excitons back from the reservoir. For the first time, we evaluate the exciton PL decay time of ~ 10 
ps. The latter process is valley-conserving and occurs on a timescale of ~ 50 ps 
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Since their isolation, atomically thin two dimensional (2D) semiconductors have attracted 
great interest due to their remarkable electronic properties and possible applications^^^. In 
particular, monolayer transition metal dichalgogenides (TMDs) show signihcant potential 
in electronic^'® and optoelectronic^'® applications. Unlike bulk TMDs, monolayers have a 
direct band gap located at two non-equivalent valleys (±i^) at the corners of a hexagonal 
Brillouin zone.^®’^^ However, perhaps the most fascinating properties of monolayer TMDs 
result from the lack of inversion symmetry and strong spin-orbit coupling. The electron 
spin and valley degrees of freedom are coupled and the resulting optical selection rules 
are valley contrasting; interband transitions in polarizations involve electron-hole pairs 
(excitons) in the ±i^ valleys, respectively and thus the valley polarization can be monitored 
through the circular dichroism.^^’^^ Polarization resolved steady-state photoluminescence 
(PL) measurements reveal a high valley polarization for monolayer MoS 2 ^^'^® or WSe 2 ^^ 
suggesting a valley relaxation time comparable to the exciton lifetime. On the other hand, 
time-resolved PL studies showed that the latter is rather short - of the order of 10 ps,^®“^^ 
setting a limit on the valley polarization lifetime. 

The observation of a robust valley polarization opens the way for applications in novel 
valleytronic devices, where the valley degree of freedom is manipulated. However, to exploit 
the potential of TMDs in this respect, a detailed understanding of the exciton PL dynamics 
and valley/spin relaxation processes is crucial. To this end, several time-resolved techniques 
have been recently applied to study the exciton dynamics in TMDs. In particular, polar¬ 
ization resolved pump-probe spectroscopy suggests that the decay of the circular dichroism 
occurs on a timescale in the range of a few to almost 100 ps.^^'^® The proposed mechanism 
underlying the loss of valley polarization is long range electron-hole exchange interaction, 
which couples the bright exciton states. 

In this paper, we investigate the exciton dynamics and inter-valley exciton relaxation in 
monolayer WSe 2 via time- and polarization-resolved PL measurements. In particular, we 
discuss the excitation density dependence of the recombination processes and of the valley 
polarization dynamics. We show that with increasing excitation power, changes in the PL 
decay are controlled by the saturation of non-radiative trap states. This is accompanied by 
a strong increase in the inter-valley scattering rate. We demonstrate that this acceleration 
affects only cold, thermalized excitons and propose that it is due to an increased exciton 
temperature or to a modihcation of valley properties of a dense exciton gas. 
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Single bulk crystals of 2H-WSe2 (the hexagonal 2H-polytype of tungsten diselenide) were 
grown using chemical vapor transport with bromine as the transport agent. Monolayer 
flakes of WSe 2 were obtained from these crystals using mechanical exfoliation. For the PL 
measurements, the sample was placed in a helium flow cryostat with optical access. Since the 
size of monolayer flakes is on the order of a few /im, the PL measurements were performed 
with excitation and collection using a microscope objective with a numerical aperture NA 
= 0.66 and magnihcation of 50 x in a back scattering conhguration. The typical diameter 
of the laser spot was approximately 1 /im. Steady-state /i-PL signal was excited with a 
640 nm laser and the spectra were recorded using a spectrometer equipped with a liquid 
nitrogen cooled CCD camera. The time-resolved /r-PL signal was excited with the frequency 
doubled tunable output of an optical parametric oscillator (OPO), synchronously pumped 
with a mode-locked Ti:Sapphire laser. The typical temporal pulse width was ~ 300 fs with 
a repetition rate of 80 MHz. The time resolved PL signal was dispersed by an imagining 
spectrometer and detected using a synchroscan streak camera with a temporal resolution of 
3 ps. Circularly polarized excitation and detection were employed to create and probe the 
non-equilibrium valley occupation. 

We start with a discussion of the steady-state PL of excitons and trions. Typical spectra 
obtained at T = 4 K from a WSe 2 monolayer as a function of excitation power are presented 
in Fig. 1(a). The peak centered at ~ 1733 meV corresponds to the recombination of the 
A-exciton, as previously reported for WSe 2 monolayers.The peak, appearing at ~ 
1710 meV is the recombination of the negatively charged exciton (trion). At lower energies 
(not shown here) a series of peaks corresponding to localized exciton states is seen.^^’^® With 
increasing excitation power, the trion emission gains in intensity with respect to the exciton 
emission. Simultaneously, the energetic separation between the exciton and trion peaks, 
AEt, increases. In Fig. 1(b), we plot AEt as a function of the intensity ratio of the trion 
and exciton emission It/Ix- We can write AEt = E^ + Ep^ where Ep and Ep are the trion 
dissociation energy (the energy required to unbind the electron from a neutral exciton) and 
Fermi energy, respectively. The latter term results from the fact that upon recombination of 
the trion, the electron occupies an available state above the Fermi level. AEp increases 
roughly linearly with excitation power, reflecting the increase of Ep, i. e., the increase of the 
number of free carriers. The plot of AEp as a function of Ip/Ix demonstrates that the ratio 
is directly proportional to the number of excess carriers. Extrapolating AE to Jt/Jx = 0 
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FIG. 1. (a) Power dependence of the cw PL spectra of a WSe 2 monolayer measnred at T = 4 K. 
Dashed lines are gnides to eye indicating increasing energy separation AEt between the exciton 
and trion PL peaks with increasing excitation power, (b) The dependence of AEt on exciton-trion 
PL intensity ratio. AEt is converted to electron concentration in the right axis. See text for details. 

yields the trion dissociation energy Ed = 19.8 ± 0.3 nieV, a value typical for monolayer 
We note that the magnitude oi Ed depends not only on the TMD material, 
but also on strain and dielectric environment.^^ The power dependence of Ep, shows that 
the number of photoexcited free carrier can be as high as 1.2 x 10^^ cm“^, for the highest 
excitation power used of 2.1 mW. 

Before turning to the time-resolved PL results, let us comment on the nature of the 
PL excitation. For steady-state experiments the excitation energy is Ep = 1937 meV, i. 
e., the excitons are created with an excess kinetic energy of E^ = Ep — Ea ~ 200 meV, 
where Ea = 1733 meV is the PL energy of the A exciton measured at low excitation power 
(see Fig. 1). The excitation is therefore below the WSe 2 band gap which is equal to about 
2.5 Although quasi-resonant (well below the continuum), the excitation results in 

the creation of a large density of free carriers. This is due to the dynamic photo-ionization 
of neutral donor states present in our n-type samples. 

In Fig. 2(a), we plot the normalized PL spectra excited with a pulsed laser at Ep = 
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1850 meV, which is close to resonance with an excited state of the A exciton.^^ In this case, 
Ek ~ 120 nieV. The presented spectra are the sum of PL signals in cr’*' and a~ detection 
polarizations. Three peaks are observed: the exciton recombination (X), the trion recombi¬ 
nation (T), and a third peak which we label as L. In Fig. 2(b), we show the dependence of 
the time-integrated PL intensity of the X peak on the excitation fluence. We find that the ex¬ 
citon intensity grows perfectly linearly over the whole fluence range used in the experiment. 
This allows us to conclude that under our experimental conditions, the exciton population 
does not exhibit any saturation from many body effects, such as biexciton formation^® or 
exciton-exciton annihilation.^^ This conclusions remains valid even at at the highest fluences 
where lattice heating induces a significant broadening and a redshift of the exciton emission 
(see Fig. 2(c)). The absence of the exciton-exciton annihilation suggests that this process 
is only efficient when assisted by exciton diffusion, which requires long exciton lifetimes and 
thus elevated temperatures.^^ 

The analysis of the PL spectra excited with OPO pulses reveals that, contrary to the 
steady state case, with increasing excitation power the trion intensity does not increase 
significantly compared to the exciton intensity, and AF^r does not increase with excitation 
power - see Fig. 2(d). This finding, corroborated on ~ 10 different ffakes, shows that in¬ 
creasing the pulsed excitation power does not result in an increased number of free electrons. 
We therefore conclude, that already at the lowest useable excitation density, the pulse power 
is high enough to ionize all the available donors. In fact, /S.Et decreases as the power is 
increased, even before the heating effects set in. We attribute this effect to a screening of 
the trion binding by the increasingly dense neutral exciton gas. The roughly linear increase 
of the trion peak with fluence (see Fig. 2(a)) reflects the creation of charged excitons from 
the carriers created by photoexcitation. 

The origin of the L peak, separated by about 54 meV from the exciton peak, was pre¬ 
viously assigned to either recombination of excitons bound on localized states^^’^®’^®, fine 
structure of the negative trion,or very recently to the recombination of biexcitons.^® As 
shown in Fig. 2(e) the L line indeed exhibits a super linear increase with respect to the 
exciton at low fluences as expected for a biexciton. However, the intensity saturates at 
high excitation powers, which together with the absence of any deviation from the linear 
dependence of the exciton intensity on excitation power, precludes the assignment of the 
L line to biexciton recombination. The fine structure splitting of the trion resulting from 
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FIG. 2. (a) Normalized PL spectra for pulsed excitation and various excitation fluences. (b) Fluence 
dependence of the integrated intensity of the exciton (X) transition compared with a calculated 
linear dependence, (c) Fluence dependence of the exciton transition energy. The redshift above 60 
/iJ/cm“^ indicates the onset of laser heating, (d) Fluence dependence of AEt- (e) The dependence 
of the integrated intensity of the L transition on the integrated intensity of the exciton transition. 
Fitted power law dependence yields an exponent of 1.3. For comparison, a linear dependence is 
plotted. 


electron-electron exchange interaction can also be excluded as the origin of the L emission. 
The recently calculated hne structure splitting of ~ 6 meV,^® is much less than the separa¬ 
tion between the T and L peaks, which is ~ 30 meV. Therefore, we assume that the L line 
originates from localized states. Further arguments in favor of this assignment are presented 
below. 

In Fig. 3(a), we compare the decay dynamics of the exciton, trion, and the L emission and 
find that the exciton decay is the fastest and that of the L line is the slowest in agreement 
with a recent report.Our temporal resolution does not allow us to resolve the rise time 
of the PL signal. All the temporal prohles can be reproduced by convoluting a Gaussian 


6 









50 100 

Time (ps) 


FIG. 3. (a) Comparison between the decay traces for the exciton, trion and the L line. IRF denotes 
the instrnment response fnnction. Its fnll width at half maximnm is 3 ps. OPO denotes temporal 
trace of the excitation pnise. Inset: the flnence dependence of the ratio between the fast and slow 
decay component (see text), (b) and (c) Temporal evolntion of the exciton and trion PL for varions 
excitation flnences. Resnlts of htting of a bi-exponential decay are plotted as lines. 


temporal profile of the OPO excitation pulse (characterized by a full width half maximum of 
about 7 ps) with a bi-exponential decay. This allows us to conclude that the thermalization 
of excitons occurs on a timescale shorter than the resolution of our measurement, i. e. 
shorter than a few ps. The fluence dependence of the PL decay traces for the exciton and 
trion transitions are shown in Fig. 3(b) and refdecays(c), respectively. The decays are htted 
with y4i exp(—f/fi) -|- ^2 exp(—f/f 2 ) and the fit is performed simultaneously to all seven 
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decay traces for fluences ranging from 4.2 to 480 /rJ/cm^. For the exciton and the trion, the 
fitting yields fi = 9 ± 1 ps and for the L transition fi = 13 ± 2 ps. The short decay time 
is related to radiative recombination and/or transfer of excitons to the non-radiative states. 
Thns, in principle ti sets the lower limit for the intrinsic radiative recombination time, in 
agreement with other reports on WSe 2 ^^ and other TMDs.^®’^° However, recent calculations 
show that the intrinsic exciton lifetime in WSe 2 at 4 K is equal to ~ 4 ps and increases 
sharply with temperature at a rate of 1 — 10 ps/K.^*^ This suggests that the initial PL decay 
could be mostly due to radiative recombination. More insight into the origin of the H time 
is provided below, in the discussion of the PL polarization. 

The fitted value of ^2 = 45 ± 5 ps of the long decay time constant is the same for all three 
transitions. We have measured several different flakes and ti always remains in the same 
range, while ^2 can be as high as 70 ps. We interpret the long decay as originating from a 
transfer of excitons from a reservoir of non-radiative (dark) defect states. Thus, t 2 represents 
the exciton lifetime in the dark reservoir. As the excitation power is increased, the reservoir 
becomes progressively filled and, as a consequence, the ratio between the amplitude of the 
short and long decays, A 1 /A 2 is increased for the exciton and the trion, as shown in the inset 
to Fig. 3(a). The observed saturation of the dark reservoir with increasing power shows that 
it is not related to states out of the radiative cone, i. e. with a high kinetic energy. For such 
states, with increasing excitation power, no saturation is expected; we should observe rather 
the opposite effect with the long decay amplitude increasing with respect to the short one. 

We are unable to resolve any variation of A 1 /A 2 for the L line. However, the decays of 
the L line could only be measured up to 60 /iJ/cm^. For higher fluences, the measurement 
was precluded by the broadened trion line and decreased intensity of the L line (see Fig. 
2(a)). 

We now turn to the discussion of the temporal dependence of the circular polarization of 
the PL signal. As shown by many authors, the circular polarization degree directly monitors 
the valley polarization of the excitons in monolayer We define the degree 

of circular polarization of the time-resolved PL signal as ^ = (/“ — -|- /^), where 

are the intensities of the PL signal in cr^ polarizations. The temporal decay of q reflects 
the decay of the valley polarization. In Fig. 4(a)-(f) we demonstrate temporal dependence 
of Q of the exciton and trion emission. It shows that during the arrival of the excitation 
pulse Q is approximately constant, with a value of about 50%. It is important to note that 
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FIG. 4. Temporal evolution of the circular polarization of the exciton PL (a)-(c) and trion PL 
(d)-(f) for various excitation fluences. (g) Valley relaxation time (tjv) extracted from the measured 
polarization decays, plotted as a function of the excitation fluence, for the three PL peaks shown 
in Fig. 2(a). 
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this initial polarization value is nearly unaffected by increased excitation fluence. g starts 
to decay after a delay of ~ 10 ps, i. e. after the excitation pulse has vanished, g decays 
roughly exponentially to a non-zero value of about 5%. This low polarization value and 
the large noise level due to weak PL signal at time delays > 40 ps make it impossible to 
conclude whether this polarization is subject to a long-lived decay^^ or rather constitutes 
an equilibrium polarization. Since its presence coincides temporally with the domination 
of the long decay of the PL signal, we suppose that this polarization results from a valley 
conserving transfer of excitons from the reservoir of dark states. Crucially, it is clearly seen 
from Fig. 4 that as the excitation fluence is increased, the decay of g is accelerated. In 
the following, we describe the currently accepted valley relaxation mechanism and within 
its framework discuss the possible origins of the observed shortening of the PL polarization 
lifetime. 

Since the excitation is quasi-resonant and cr’''-polarized, due to the valley contrasting se¬ 
lection rules, the excitons are formed uniquely in the +K valley.However, since the ex¬ 
citons are created with a non-zero center-of-mass momentum k resulting from the excess en¬ 
ergy Ek, an efficient depolarization process takes place. As proposed by many authors, 
the underlying mechanism is related to the long-range part of the electron-hole exchange 
interaction, which can be treated as an effective magnetic held f2 around which the valley 
index precesses. Momentum scattering events lead to a randomization of the held direction 
and thus the valley relaxation by the mechanism proposed for quantum well excitons.We 
assume a strong scattering regime, in which flTp 1, where Tp is the momentum scattering 
time. Under such conditions the intervalley scattering time is where the 

brackets denotes averaging over the distribution of exciton energies.Since fl ~ 
the valley polarization, and thus the circular polarization of PL, exhibit a strong dependence 
on the excitation energy.In particular, for the values of Ek from our experiments, an 
initial depolarization time on the order of 10 fs is expected.^® After thermalization and sub¬ 
sequent cooling, which occur during hrst ~ 0.5 ps after excitation,the average k value of 
the exciton population is reduced and, consequently, the valley depolarization slows down. 

In order to reproduce the measured temporal dependence of g and to extract the exper¬ 
imental values of we employ a simple rate equation model for the exciton populations 
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in the ±K valleys: 
dN^ _ N+-N- 

dt tree ^ tiv 

where tree is the recombination time. Taking into account two exciton populations, neglects 
the fact that the exciton exhibits a richer fine structure with two additional dark configu¬ 
rations with parallel electron and hole spins. However, since the long-range electron-hole 
exchange couples only bright exciton states,only these two states are involved in the 
inter-valley scattering process, which allows us to neglect the dark states. In the case of 
the negatively charged trion, there are six bright spin/valley configurations split by both 
electron-hole and electron-electron exchange interactions.^^ Since the splitting between the 
ground and the first excited trion configuration is about 6 meV,^® which is an order of mag¬ 
nitude larger than the thermal energy at the temperature of the experiments, we assume 
that only the ground configuration is populated and neglect the higher energy ones. 

To account for the ultra fast initial depolarization, we shift the time origin by about 
10 ps, i. e., to the time delay when the exciting OPO pulse has vanished and p starts to 
decay. Furthermore, we choose the initial conditions as iV^(O) = p and N~{0) = 1—p, where 
2p — 1 is the initial value of p. The solution to Eq. 1 yields for the temporal dependence 
of the polarization; g{t) = {2p — 1) exp(—2t/tiv). Note that is twice the polarization 
decay time. The experimental data is fitted with this formula plus a constant background 
evaluated as a mean value of g for decays between ~ 75 and 85 ps, accounting for the 
background polarization. This procedure allows to evaluate p and We find that for the 
exciton p = 0.72 ± 0.07, approximately constant over the whole range of excitation fluences. 
For the trion, p decreases monotonically from 0.75 to 0.65 as the fluence is increased. The 
fluence dependence of fiv plotted in Fig. 4(g) reveals that with increasing fluence, intervalley 
scattering time decreases from ~ 50 ps to ~ 20 ps for both the exciton and trion. 

A natural question arises whether the acceleration of the intervalley scattering can be 
attributed to the heating of the lattice by the laser pulse. Increased scattering via an in¬ 
creased acoustic phonon population is ruled out since it would require heating of the lattice 
to temperatures higher than 100 At the temperature of the experiments scattering on 
ionized impurities is expected to be the dominant momentum relaxation mechanism. An 
increased lattice temperature would also lead to a reduction of the band gap and, conse¬ 
quently, a larger since we keep the excitation energy constant. This would lead to a 
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decrease of the initial polarization via an increase of both f2 and Tp. However, as seen from 
the temporal traces of g presented in Figs. 4(a)-(f), polarization valnes at short time decays 
are approximately independent of the excitation fluence (a fact corroborated by a very weak 
dependence of the htted p as a function of fluence). Moreover, the dependence of the exciton 
transition energy on fluence (see Fig. 2(c)) demonstrates the absence of any lattice heating 
at least up to about 60 /rJ/cm^. 

Another mechanism leading to the acceleration of intervalley scattering could be related 
to a screening of the Coulomb potential of the ionized impurities. As shown for intrinsic 
GaAs quantum wells, this mechanism was found to be responsible for the decrease of spin 
relaxation time with increasing excitation density as a result of an increase in However, 
this mechanism would also affect the initial polarization and therefore can be ruled out in 
our case. 

Our experimental results demonstrate that the increased intervalley scattering rate con¬ 
cerns only cold excitons, i. e., with a small center-of-mass momentum. We propose that 
the effect leading to the increased valley relaxation rate can be related to the heating of the 
exciton system. Indeed, larger average kinetic energy results in a larger f2. The heating 
could occur as a result of an increased exciton-exciton scattering. Elevated quasi-particle 
temperature explains also the small decrease of the initial trion polarization; increased trion 
temperature leads to a small occupation of the excited trion conhgurations^^’^® thus opening 
additional intervalley scattering channels. As shown by the temperature dependent time- 
resolved Kerr rotation studies, enhanced intervalley relaxation requires exciton tempera¬ 
tures higher than 30 We remark that at such high temperatures we expect a signihcant 
increase of the exciton radiative lifetime, in agreement with a recent theoretical work.^° 
However, the measured ti does not exhibit any dependence on fluence suggesting that H 
contains a signihcant contribution from non-radiative processes. 

Another possible origin of the increased intervalley scattering rate can be attributed to the 
presence of a dense exciton gas. We calculate the huence dependent average exciton-exciton 
distance by assuming that about 1% of the photons are absorbed by the WSe 2 monolayer^® 
and that each absorbed photon is converted into an exciton. This crude assumption gives 
a lower (upper) bound for the average exciton-exciton distance (exciton density). The cal¬ 
culated values are given in the upper axis of Fig. 4(g). It shows that at the highest huence 
the average exciton-exciton distance becomes comparable to the exciton Bohr radius of 1.3 
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Under such conditions, screening and/or many body effects are expected to inffuence 
the exciton valley dynamics. 

In conclusion, we have studied the exciton and trion PL in monolayer WSe 2 . The intensity 
of the exciton PL shows no signs of saturation up to the highest excitation densities where 
heating effects are already present. This demonstrates that non-linear processes such as 
exciton-exciton annihilation or biexciton formation are inefficient at low temperatures and 
for quasi-resonant excitation. In addition, we have shown that dark non-radiative trap states 
play an important role in determining the PL decay dynamics. However, at high excitation 
densities these states become saturated. Finally, the measured intervalley scattering rates 
of both excitons and trions are significantly increased with increasing excitation density 
which may have negative effects upon the operation of valleytronic devices. Importantly, 
the increased intervalley scattering rate occurs only for cold excitons and does not inffuence 
hot excitons - created with a high kinetic energy. We attribute this effect to either heating 
of the exciton system by increased exciton-exciton scattering or to the presence of a dense 
exciton gas leading to a modification of the exciton valley properties. 
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